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Abstract 

The magnetoresistance oscillations in GaAs/Al^-Gai-ajAs heterostructures induced by millimeter- 
wave radiation recently observed by Zudov et al. [Phys. Rev. B 64, 201311 (2001)] is theoretically 
reproduced by introducing a model based on the finite temperature Fermi liquid theory of dc con- 
ductivity and the electron reservoir hypothesis. The obtained oscillation patterns show excellent 
agreement with the experimental result by Zudov et al. and are independent of the polarization 
of the radiation field in accordance with the experimental observation by Smet et al. [Phys. Rev. 
Lett. 95, 116804 (2005)]. 
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The new class of magnetoresistance oscillations in a high-mobility two-dimensional elec- 
tron system (2DES) in a GaAs/Al^Gai_a;As heterostructure subjected to weak magnetic 

21] yl and by Mani et al. 



:ields and millimeterwave radiation found by Zudov et al. [1 
has revealed again the fascinating nature of the 2DES 



111 ], evok- 



ing the discoveries of the quantum Hall effects (QHE) by von Klitzing 12| and by Tsui 
et al. [iS^. The period of these new oscillations is governed by the ratio of the mil- 
limeterwave to cyclotron frequencies, and the minima of the oscillations are character- 
ized by an exponentially vanishing diagonal resistance. When the millimeterwave radi- 
ation is turned off, the magnetoresistance shows the well-established SdHvA oscillation 
whose period is governed by the ratio of the chemical potential to the cyclotron fre- 
quency. Hence this magnetoresistance oscillation is apparently induced by the illumination 
with millimeterwaves. Although various different theoretical models have been proposed 
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3J], one crucial 



question remains to be solved. The experiment by Smet et al. shows that the resistance 
oscillations are notably immune to the polarization of the radiation field. This observation 
is discrepant with these theories and seems to cast doubt on the validity of the theoretical 



models so far proposed 



361 ]. In these circumstances it seems necessary to explore other 



theoretical possibilities to explain the phenomena. 

Meanwhile, an extremely interesting observation was made by Holland et al., who mea- 
sured the long-wavelength magnetoplasmon dispersion in a high mobility 2DES realized in 
a GaAs/Al2.Gai„^As heterostructure quantum well 37|]. Using the coupling between the 
plasmon and THz radiation, they explored a wide range of filling factors for a fixed value of 
the wave-number vector to obtain an explicit filling-factor dependence of the dispersion. The 
observed dispersion shows a clear filling-factor dependent plateau-type dispersion, violently 
deviating from the well-established semiclassical formula at first glance. Their discovery 
suggests a previously unknown relation between the magnetoplasmon and the integer QHE 
(IQHE). The full theoretical explanation of the phenomenon was given very recently 38] 



by adopting the electron reservoir 



decades ago to explain the IQHE 



hy 



39, 



3othesis (ERH) that was introduced more than two 
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I4I]. That the experimental discovery by Holland 
et al. can be explained by adopting the ERH strongly indicates the significance of taking 
into account the electron reservoir in order to investigate the behavior of the electrons in the 



GaAs/A^Gai_2:As heterostructure 2DES. In this paper a Fermi liquid model based on the 
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ERH is proposed to explain the new class of the magnetoresistance oscillations within the 
framework of the standard quantum theory of electrical conductivity. The derived formula 
gives the magnetoresistance oscillation patterns which show excellent agreement with the 
experimental data observed by Zudov et al., including the SdHvA oscillation part. 

The Fermi liquid theory of electrical conductivity was originally formulated by Eliash- 
berg. The core of the theory is the analytic continuation of the finite temperature current 
correlation function with respect to the Matsubara frequency to obtain the retarded real- 
time current resp onse function, which directly yields the conductivity. Here the formulation 
given in Ref.j42l] is applied to an electron gas confined in the xy-p\&ne and subjected to a 
magnetic field B = (0, 0, B) and a millimeter wave whose angular frequency is u. The polar- 
ization of the radiation field is arbitrary. The electron charge and effective mass are denoted 
by — e and m, respectively. The interactions between the electrons and impurities or de- 
fects are treated perturbatively by the Bethe-Salpeter (BS) equation for the vertex function. 
The quasi-electron spectrum is assumed to be e = tujJc{M + 1/2) + {g* hbB /2)sgYia = Sua-, 
where M is the principal quantum number for the Landau levels, Uc = eB/mc is the cy- 
clotron frequency, g* is the effective (^-factor, fiB = eh/Airmoc is the Bohr magneton with the 
electron rest mass mo, and the function sgn gives the sign of the spin variable a. The Lan- 
dau levels are degenerate and there is an additional quantum number p, which corresponds 
to the a;-component of the electron momentum. The allowed range of this momentum is 
bl < Pmax = neB/hc. 

Although a microscopic model for the scattering mechanism responsible for the conductiv- 
ity may be introduced via the proper vertex function in the BS equation, here we do not use 
a specific model but simply assume that the scattering mechanism is independent of M and 



a. Then applying the theory formulated in Ref. 
the general expression of conductivity can be found as 



42| to the two-dimensional quasi-electrons. 



4m^ 

a M=0 



E E ^ / ^^ dPP'^Gf^aiP^u)GUp,^) 

a M = J -Pmax 

X |l + lReAn(p,a;)|, (1) 

where f{uj) = (1 + eyi\i{(3hjj))^^ , (G^) is the retarded (advanced) Green function, and 
An is the vertex function. The aim of this paper is to examine the general consequence of 
this conductivity formula when the ERH is adopted. 
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The fundamental assumption of the Fermi hquid model, i.e., the Fermi liquid hypothesis 
(FLH), can be mathematicall y fo rmulated in terms of the retarded and advanced Green 
functions, and G^. In Ref. 42| it has been shown that by virtue of the FLH the product 
G^G^ in the above conductivity formula can be written as 

,2 



rra 



7 



6 (^uj — h ^(e — /i)) + {non-singular term} , 



(2) 



where a is the wave function renormalization constant, and 7 is the damping parameter, 
which should satisfy the condition phj -C 1, and At is the chemical potentia. 



In the theories so 
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34l |. the effects of the millimeterwave radiation only on the 2DES 



were considered. However, if there is an electron reservoir, the effects of the radiation on the 
electrons in the reservoir should also be taken into account. At present the mechanism of the 
electron reservoir is not known. Here electrons in the reservoir are assumed to be in bound 
states with binding energy —Eres- Since the amount of the energy that an excited electron 
can receive from millimeterwave radiation is hu, the condition that the electron can join the 
2DES should be f}Wc/2 + E^es < which can also be written as 5 < {2mc/ he){hi' — E^es) = 
Be- As the chemical potential is the minimum free energy to add an electron to the system, 
the emergence of such a process may be described by introducing another singularity in the 
retarded Green function at hu. This singularity of the Green function may be expressed as 
an effective chemical potential in the expression such as (2). These considerations lead to 
the following expression for G^G^: 

,2 



G^{u!)G^{u;) = — -6 {ljJ — h ^{smol — Vi)) + {non-singular term} 



1=1,2 

when the radiation is on, and 
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(3) 



72 



6 (lj — h '^{sMa — V2)) + {non-singular term} 



(4) 



when the radiation is off. Here we have defined rji = hu, 772 = /i, and regularized step 
functions = (1 + expr(i? — Bc))^^ and = (1 + expr(i?c — B))^^. The limit r — 00 
corresponds to the sharp cut-off. By virtue of the delta functions the integration over uj in 
(1) for the two terms in (3) becomes trivial. Since the damping takes place for non-zero 
value of the momentum, the main contribution to the p-integration comes from the vicinity 
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of p = 0. That is, we can safely approximate 7^ ~ 7oiP^- Then the conductivity formula (1) 
yields 

^ Awiv? ^ C/T ^ {^i/c^(^Mq - hv)Oi + \2fu{eMa - 1^)02 + Ao} (5) 
a M=0 

when the radiation is on, and 

-e^h ^ eB 



5Z ^^2/c.(eMa -/i) + Ao} (6) 



Aum? ^ ch 



when the radiation is off. The quantities Aj = a^(87r7oj)^^ {1 + h^^KeAu} for i=l and 2 
depend on the microscopic mechanism of the scattering process responsible for the conduc- 
tivity. The contribution from the non-singular term is denoted by Aq. For simplicity, we 
assume Ai = A2 = A. Using the explicit form of /(cu), the conductivity can be expressed as 

a.. = — {W,e^ + W^e^ + iB} = —Won (7) 
m m 

when the radiation is on, and 

= —{W2+ m = —WoFF (8) 

mm 

when the radiation is off. Here ^ = {—e/Air'^mcX) ^j^^ Aq is assumed to be a constant, and 
Wi's are given as 

n _|_ pP{£Ma-V- 

a M=0 

This conductivity formula is derived within the theoretical framework of the linear response 
approximation, in which the electric field E = [E^, Ey) is an externally controlled small per- 
turbation. In the measurement by Zudov et al. the current Ix is measured by controlling 
the electric field E^-, while the current ly as well as the external electric field Ey are kept 
zero. Therefore, the resistivity R^x observed in their measurement should correspond to 
1 / (^xx in this theory. The significance of considering the boundary conditions on the current 



' Z-^ Z-^ |1 -|- e/3(eMa-»?i) 1 |i -|- e-/3(£Mc«-'?i)| ^ ' 



expectation values has been discussed in Ref. 



4l| . where the IQHE formula is derived from 



a microscopic many-body Hamiltonian without depending on the linear response approxi- 
mation. Thus, within the theoretical framework of the present formulation the resistivity 
corresponding to Rxx in Ref. [ll is given as 



TTl 

Rxx = (J~x = ^^ON or OFf(^, T, fX, u) . (10) 



The function Wqn depends on three adjustable parameters r, and ^, while Wqff depends 
only on ^. These parameters will be adjusted to optimize the fit of the theoretical curve 
with the experimental data. 

The chemical potential of 2DES without millimeterwaves is given as 



, (11) 

m 

where n is the electron number density in the limit of very weak magnetic field. This gives a 
good approximation for the range of the parameters in the experiments of Ref . [l| . It should 
be noted here that the above approximation gives the classical Hall effect instead of the 
IQHE. This explains the experimental observation reported in Ref.jsj. 

The 5-dependence of Wq^ for 0.2 < B < 2 kG is shown by a solid line in Fig. 1, using 
the values of the physical parameters given in Ref. i.e., T = 0.5 K, n = 2 x lO^^cm"^, m = 
0.068mo, and / = 100 GHz. The adjustable parameters are chosen such that B^ = 3.0 kG, 
r = 0.01, and ^ = 1.8 x 10~^. These values indicate Eres — 0.16(meV). The experimentally 
measured R^x taken from Fig. 4 in Ref. Ij is also shown by a dotted line for the same B 
values. The oscillation pattern of Wq^ shows excellent agreement with the experimental 
result. In Fig. 2 the -B-dependence of Wq^ is shown by a solid line for 0.2 < B < 3.6 kG, 
using the same values for the parameters. The experimentally measured R^x taken from 
Fig. 4 in Ref. |l| is also shown by a dotted line for the same B values. The theoretical 
pattern shows excellent agreement with the expermental curve. The phase of the oscillation 
observed in the experimental curve is shifted toward positive B direction. The shift becomes 
prominent for B > 2.2 kG. Theoretically, an additional term to the quasi-electron energy 
spectrum, i.e., the proper self-energy, may cause such a shift. The magnetoresistivity data 
given in Fig. 1 of Ref. [l] show a peak due to the magnetoplasmon resonance around B ~ 2.4 
kG. The existence of magnetoplasmon resonance peak in the magnetoresistivity in Hall bar 
structures was first observed by Vasiliadou et al. [43| in the microwave photoconductivity 
measurement. The magnetoplasmon resonance gives considerable contribution to the two- 
particle Green function, which is directly related to the proper self-energy of the single 
particle excitation spectrum via the finite temperature Ward-Takahashi relation 
Fig. 3, the function W^JpF plotted for 0.2 < B < 3.6 GHz. It exhibits clear SdHvA 
oscillation observed experimentally. 

In conclusion, the 5-dependence of the oscillatory patterns of the millimeterwave induced 
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35|. Since the model 



magnetoresistance oscillations observed by Zudov et al. is almost perfectly reproduced 
from our theoretical model based on the FLH and the ERH. Furthermore, it should be 
noted that the present model is independent of the polarization of the raditation field in 
perfect accordance with the experimental observation by Smet et al. 
does not depend on any specific mechanisms of the electron scatterings that cause electric 
resistivity, this new magnetoresistivity oscillation seems to be a universal consequence of the 
Fermi liquid nature of the 2DES and the existence of an electron reservoir in GaAs/AlGaAs 
heterostructures. 

In view of the fact that the present theory can reproduce excellently the experimental 
data including the immunity to the polarization of the radiation field, the next task should 
be to reveal the microscopic mechanism that realizes the electron reservoir. At present there 
seem to be no direct experimental clues to that question. Nevertheless, it should be noted 
that there has been no direct experimental evidence against the existence of an electron 
reservoir either. In other words there has been no direct experimental evidence that the 
electron number is fixed. 

As the essential theoretical conjecture of this work, i.e., Eqs. (3) and (4), can be used 
to calculate various physically measurable quantities, the validity of the model presented in 
this work may be further tested. 
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Fgiure captions 



Fig. 1 

Tlie 5-dependence of tlie function Wq^ and tlie experimentally measured R^x for 
0.1 < B < 2 when the system is illuminated with radiation. The parameters are: / =100 
GHz, T = 0.5 K, n = 2 X lO^^cm'^ . (a) Theoretical W^^ defined by Eq. (7) is plotted by 
the solid line, (b) Experimentally measured R^x taken from Fig. 4 in Ref. l| is plotted by 
the dotted line. 



Fig. 2 

The -B-dependence of the function Wq^ and the experimentally measured Rxx for 
0.2 < B < 3.6 when the system is illuminated with radiation. The parameters are same 
as Fig. 1. (a) Theoretical Wq^ defined by Eq. (7) is plotted by the solid line, (b) 
Experimentally measured Rxx taken from Fig. 4 in Ref. [l| is plotted by the dotted line. 



Fig. 3 

The -B-dependence of the function H^oFp for 0.2 < B < 3.6 when the system is not 
illuminated with radiation. It shows a typical SdHvA oscillation. 



11 




0.5 1 1.5 2 



B{kG) 



Fig. 1 




Fig. 2 



2 



1.5 



0.5 








B(kG) 



Fig. 3 



